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Abstract. Two effective methods for writing the dynamical equations for non-holonomic 
systems are illustrated. They are based on the two types of representation of the constraints: 
by parametric equations or by implicit equations. They can be applied to linear as well as 
to non- linear constraints. Only the basic notions of vector calculus on Euclidean 3-space 
and on tangent bundles are needed. Elementary examples arc illustrated. 
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1 Preamble 

The classical theory of non-holonomic dynamical systems, even in recent times, is treated in 
a growing number of papers. Most of them require the use of modern differentiable and algebraic 
structures which are not familiar to non- mathematicians working on concrete applications. On 
the other hand, several papers are dedicated to the analysis of special non-holonomic mechanical 
systems, quite interesting but treated with ad hoc methods. These are the reasons why I think 
useful to propose a ready to use approach to the dynamics of non-holonomic systems, requiring 
the knowledge of the basic notions of the vector calculus on the Euclidean three-space and on 
tangent bundles only, and avoiding the use of cotangent bundles (Hamiltonian framework) and 
jet-bundles. 

In the present paper I illustrate two different, general, effective and concise methods for 
writing the dynamical equations of a given non-holonomic system. These methods correspond 
to the two ways of representing kinematical constraints: by parametric equations or by implicit 
equations. They can be applied to linear as well as to non-linear constraints and lead to different 
(but obviously equivalent) dynamical systems. 

In order to make this paper self-contained, a straightforward approach to the Gauss principle 
and to the Gibbs-Appell equations is illustrated. Our starting point will be the well under- 
standable Newton dynamical equations for a system of massive points. 

The tutorial character of this paper does not exclude the presence of some novelties. 

Many articles and books have been consulted in writing this paper. To cite all them would 
make the list of references quite long. On the other hand, such a long list would not fit with 
the purposes of this paper. Anyway, I must mention the excellent recent books |2j [3], from 
which I have got many useful hints. They are very well readable and advisable to non-experts 
which would like to go into this matter in more depth, and to know its present developments 
and applications. 

*This paper is a contribution to the Vadim Kuznetsov Memorial Issue "Integrable Systems and Related Topics". 
The full collection is available at |http://www.emis.de/journals/SIGMA/kuznetsov.html| 
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2 Introduction 

Let us consider a mechanical system with a well defined configuration manifold Q n . The di- 
mension n of Q is the number of degrees of freedom of the system. Let us denote by q = (q l ) 
any (Lagrangian) coordinate system on Q and by (q, q) = (q l ,q l ) the associated coordinates on 
the tangent bundle TQ. This tangent bundle is the space of the kinematical states. 

A kinematical constraint is given by a subset C C TQ. A special case is that of regular 
constraint: 

Definition 2.1. A kinematical constraint C C TQ is said to be regular if 

• C is a submanifold of dimension n + m, m < n; 

• for all q £ Q, F q = C H T q Q is a submanifold of dimension m; 

• the restriction to C of the tangent fibration tq : TQ — > Q is a surjective submersior^. 
The constraint is linear if each is a subspace (see Fig. 1). 



Space of kinematical states TQ 








Constraint C 
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Configuration manifold Q 
Figure 1. A regular constraint is a fibration onto Q. 

A regular constraint can be represented in two ways: 

• Parametric representation: it is described by m equations, 

<f = V%,*), (2.1) 

where z = (z a ), a = 1, . . . ,m < n are called 'parameters. Note that (q, z) can be interpreted 
as local coordinates on C. 

• Implicit representation: it is described by r = n — m independent equations^], 

C a (q,q) =0, a=l,...,r. (2.2) 

x If this notion is not understood by the reader, he can look at the equivalent conditions (|2.3jl and (|2.4|l below. 
The definition of regularity is taken from [10] and [6]. In an extension of this definition is given: it requires the 
existence of a submanifold Qi C Q such that C is a submanifold of TQi and the restriction to C of the tangent 
fibration is a submersion. For our purposes we do not need to consider this more general case. 

2 This means that the differentials dC a are linearly independent at each point of C. 
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In these two representations, the regularity of the constraint is represented by conditions 
rank [ip l a } nxm = m, rank [Cf] nXr = r, (2.3) 



respectively, where 

tf'-M C^ 8 -^ (24) 



Note that the regularity conditions may be not fulfilled at certain states, that we call sin- 
gular states. They always occurs, for instance, for non-linear homogeneous constraints (see 
Remark 17. 5p . 

In the following, with the exception of special remarks, the constraints are assumed to be non- 
linear. However, in the case of linear constraints, we shall assume - without loss of generality - 
that the functions t^ l (q,z) are linear in the parameters, so that ip l = ^p t a (q)z a , or that the 
functions C a (q, q) are linear in the q, so that C a = Cf{q) q l . 

The leading idea of the first method is to consider the parametric equations (12. ip of the 
constraint as a first set of dynamical equations, to be completed by a second set of equations of 
the kind z a = Z a (q,z): 

q l = ^(q,z), 

z a = Z a (q,z). (2.5) 



This is a first-order system of ODE's. Once the initial conditions (go, zq) are fixed, they give 
a unique actual motion of the mechanical system. By actual motion we mean a physical 
motion obeying to the constraints. 

The explicit expressions of the functions Z a (q,z) depend of course on the given functions 
ifj l (q, z). But they depend also on the dynamical principles we assume. 

The Newton dynamical equation for a material point will be the only physical principle on 
which we shall base our approach. Furthermore, we shall assume that the constraints are ideal 
or perfect, according to a suitable mathematical definition: this means to accept a certain 
constitutive condition on the constraint as a postulate. In this way, the Gauss 'principle 1 will 
follow as a ' theorem'' from the Newton equations of a system of particles. 

The second method presented here is based on the implicit representation of the con- 
straints. It is a revisitation, with improvements and simplifications, of the well-known Lagrange- 
multipliers method. In this context we get the explicit expression of the reactive forces. This 
is useful, in the concrete applications, for measuring the stress that the constraints have to 
support. 

These two methods lead to different first-order dynamical equations i.e., to different vector 
fields: 



• The first method (parametric representation) produces a vector field Z on the constraint 
manifold C. The integral curves of Z give all actual motions - see Fig. 2. 

• The second method (implicit representation) produces a vector field D on the whole TQ 
but tangent to C. Only its restriction to C has in fact a physical meaning i.e., only its 
integral curves which start from a point of C (and which will lie on C) represent actual 
motions - see Fig. 3. As a byproduct, this method gives the explicit expressions of the 
reactive forces, which can be estimated along any actual motion. 
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TQ:(q,q) 




Q--(q) 

Figure 2. Parametric representation of C - The 
vector field Z. 

3 Ideal constraints 



Q-Ui) 

Figure 3. Implicit representation of C - The 
vector field D. 



At a microscopical level, a mechanical system is made of a collection of material points {P u ,m u ), 
v = 1, ... ,7V. The position vector r u of each point P v in the Euclidean three-space is a func- 
tion r y {q) of the chosen Lagrangian coordinates. A motion of the system is then represented 
by a time-parametrized curve q l = q l {t) on Q, and at each configuration q £ Q all possible 
velocities of the points are given by 



dq 



T 1 



(3.1) 



A kinematical state of the system is the collection of all pairs position- velocity (r^, v^) of the 
points. The collection of all possible states is then the tangent bundle TQ of Q. At any fixed 
state the accelerations & u are given by 



d 2 r u . { ■ div dif_ 
dq l dq^ dq' 1 dt 



(3.2) 



Let us consider the parametric representation (|2.1j) of the constraint. Then the velocity and the 
acceleration of each point, compatible with the kinematical constraint, are 



dq 1 



Let us write these equations in the more compact form, according to the above-given notation: 
v v (q,z) = div v ip\ 

8i u (q, z, z) = dijT v W + div v {i>) V + z a ) . (3.3) 



As we shall see, the additional parameters z play a crucial role. Note that, at any fixed 
configuration q £ Q, the parameters z determine all kinematical states (q,v) compatible with 
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the constraints. Intuitively, the parameters z determine a small displacement from a state (q, q) 
to another close state (q,q') (with the same configuration q). For a precise definition, let us 
decompose the acceleration (|3.3p into the sum 



a„ = a Qu + a au z a , (3.4) 

where 

&0is = dijr u VV J + d { x v 4>j V> J , a au = diV u ip l a . (3.5) 

We observe that the vectors aoj, depend only on the state of the system. We say that 
Definition 3.1. The second vector in (|3.4p 

w v = a av z a , (3.6) 
represents the virtual displacements of the point P u at the given kinematical state. 

Remark 3.1. Note that the parameters z a span all possible virtual displacements at a given 
state. The mechanical meaning of virtual displacement, which is strictly related to that of 
ideal constraint (Definition I3.2|) would require a detailed discussion. We can skip it simply by 
accepting equation (|3.6f) as a mathematical definition, since it will be justified first, by the fact 
that with such a definition the Gauss principle becomes a consequence of the Newton dynamical 
equations (Theorem I4.ip and second, by the fact that the reactive forces of ideal constraints are 
not dissipative (Remark 13. 3p . 

It is customary to associate the intuitive idea of 'virtual displacement' with that of 'virtual 
velocity', as a limit of a 'small' displacement between two configurations' of the system. Instead, 
within the present context, a 'virtual displacement' is a 'small' displacement between kinema- 
tical states (configurations plus velocities), so it is associated with the intuitive idea of 'virtual 
acceleration'. This viewpoint is in fact coherent with the philosophy of the Gauss principle, 
which deals with accelerations. 

We assume for the dynamics of each point (Pj,,mj,) the Newton equation 

m u a v = A y + Rj,, (3.7) 

where A„ is the active force (due to external fields and internal interactions) and Rj, is the 
reactive force: it has the role of making the constraint fulfilled. 

Remark 3.2. The idea of 'reactive force' arises from the Newtonian philosophy, according 
which any action deviating a point from the uniform rectilinear motion (in an inertial reference 
frame) is a 'force', mathematically represented by a vector. Thus, the presence of a kinematical 
constraint must be represented by a vector, called 'reactive force', to be summed to the 'active 
force', which in turns represents the action of fields present in the space and independent from 
the constraints (gravitational, electromagnetical, centrifugal, Coriolis, etc.). 

Definition 3.2. Non-holonomic constraints are said to be ideal or perfect if 

R„ • w„ = (3.8) 



for all virtual displacements w u . 
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We consider equation (|3.8p as a constitutive condition for the constraint: it says which kind 
of reactive forces the constraint is able to supply in order to be satisfied along any motion. It is 
straightforward to check that for linear constraints, or simply for holonomic constraints (which 
do not involve velocities), equation (|3.8p reduces to the classical virtual work principle. The 
validity of such a constitutive condition is a matter of theoretical and experimental analysis of 
the behavior of the mechanical system one is dealing with. 

Theorem 3.1. Let 

V 

be the Lagrangian reactive forces. Then the definition (13. 8p of ideal constraint is equivalent 
to the following equations, 

Ri C = 0, (3.10) 
Ri = \ a C?, (3.11) 

in the parametric and in the implicit representation, respectively. 

Equation (|3.1ip means that the components Ri are linear combinations of the functions 

Proof. Put equation (|3.6p and (|3.5|) in equation (|3.8j) . 
o = • *<*> z a = ]T • div v ipi z a = Ri < z a . 

V V 

This proves equation (|3.10p . By differentiating the identity C a (q, ijj(q, z)) = 0, we get the 
following relations between the two representations, 

CfC = 0, diC a + C^i = 0, (3.12) 

where 

fina P) 

Ct = W , ^=8^, dl = g-. (3.13) 

Due to the first identity (|3.12p . equation (|3.10p is then equivalent to (|3.1ip ■ 

Remark 3.3. Equations (|3.13P i and (|3.1ip show that the ideal constraints do not dissipate ener- 
gy: the power of the reactive forces is zero (we are dealing with time-independent constraints). 
See 0. 

Remark 3.4. By a well-known process, we pass from the microscopical level to the l macrosco- 
picaV one i.e., to the Lagrange equations 

d fdK\ 8K . t . 

where 

K = g 9ij Q t( l j , 9ij = ^Z m v diY v • djY v . (3.15) 

V 

At the right hand side of the Lagrange equations we find the sum of the active Lagrangian 
forces 

V 

and the reactive Lagrangian forces (|3.9p . 
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4 The Gauss principle 



At the microscopical level we introduce the quantity 

^ 2 

2 ^ " V TTl^ 

The active forces A„ are known functions of the state (g, q). Thus, due to the parametric 
equations of the constraints, G becomes a function of (q,z). Moreover, even for active forces 
depending on the velocities, A„ does not depend on i, 

|^ = o. 

Thus, due to equations (|3.4p and (|3.5p . along any motion satisfying the constraints we have 

%f = £ ™» - • §S = E m * - ^) • a - (42) 

Furthermore, 

d 2 G 5a 



^ m„ ^-r^ • a Ql/ = ^ m„ • a ai/ = ^ m v • djr u = g i:j ip % Q ip j p . 



dz a dzP ^ v dz? 

V V V 

Then, if we introduce the functions 

G a p = gtj </i ^j, (4.3) 
we get 

d 2 G 



dz a dzP 



G aP . (4.4) 



Since the matrix [ijj l a ] has maximal rank, the symmetric matrix [G a p\ is regular and positive- 
definite as well as [gij]- 

Theorem 4.1. Assume the Newton equations m u a. p = A u + for each point P v . Then, at 
any state of any actual motion the quantity G takes a minimal value (Gauss principle) if and 
only if the constraints are ideal. 

Proof. Write the Newton equations in the form 
m v I a„ - — ) = R 



m 

Then, due to equations (|3.4p . (|3.5p and ([37 

V V 

This shows that for ideal constraints, see equation (13.10p . the Newton equations imply 
dG 



dz a 



0, (4.6) 



at any state along any actual motion. Due to equation (|4.4p . being [G a p\ positive, at the statio- 
nary states for which equation (|4.6p holds, the function G has a strong minimum, (ii) Conversely, 
assume that the Gauss principle holds true. Then equation (|4.6|) is satisfied, so that from 
we get Ri ip l a = 0. This means that the constraint is ideal (Theorem 13. ip . 
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Remark 4.1. The vector A u /m u is the acceleration of the point P v in a free motion, free 
from the constraints. Let us denote it by a u . As a consequence, the function G can be also 
defined as 

G = ^M*»-4) 2 , (4-7) 

V 

and Theorem 14.11 can be reformulated as follows: 

Theorem 4.2. Let r u {t) and rt{t) be two motions of the system P v such that for t = to the 
corresponding states coincide i.e., 

r v (t ) = rl(t ), Vv(to) = v£(to). 

Assume that r[(t) is a free motion. Then, at this state, and for any motion compatible with 
ideal constraints, the actual accelerations a. u (to), are such that G takes a minimal value. 

Remark 4.2. For any arbitrary motion, 

dr v mi dr u H ■ dr u , 4 

v v = tt-£ Q , a u = qq J + . 

oq l oq l oq^ aq l 

Then, at any fixed state compatible with the constraints we have 

^-4 = Of -<?})■ 

Due to the definition (|3.15p of gij, from (|4.7p we get the so-called Lipschitz expression of the 
function G: 

G=\ gij (4.8) 
Note that in this expression the Christoffel symbols are not involved. 

5 The Gibbs— Appell equations 

Let us go back to the definition (|4.1h of the function G. If we introduce the functions 

5 = i^ m "^ ^ = ^E — A - S 2 = J2^-^, (5.1) 

v v v 

then we have the decomposition 
G = S + Si — S2. 

The function S is called the energy of the accelerations. We observe that 

and that, due to the second equation (|3.3p and the definition (|3.16p of active Lagrangian force, 

V 



Thus, 

dG OS 



dz a dz c 

Due to the Gauss principle (Theorem 14. ip . this proves 
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Theorem 5.1. The Gauss principle is equivalent to equations 
81 

W = A.. <M> 
where the function 

S(q, z,z) = 2^2 m " a l ( 5 - 3 ) 

V 

is determined by the expression of the accelerations, and 

A a = A^ a . (5.4) 

Remark 5.1. Equations (|5.2|) are the celebrated Gibbs Appell equations. The equivalence 
between these equations and the Gauss principle is highlighted within the framework of the 
parametric representation (|2.1|) of the constraints. 

Remark 5.2. The quantities A a can be computed by writing the the virtual power of the 

active forces: 

w = A - • w - = A » • dir » < * a = Ai ^ iQ = A « * a - ( 5 - 5 ) 

V V 

6 The explicit form of the Gibbs— Appell equations 

Both sides of the Gibbs-Appell equations (|5.2p are functions of (q, z, z). Let us solve them w.r.to 
the variables z a . To this end, it is crucial to observe that by using the second equations (|3.3[) 
we get for the function S fj5.3j) the expression 



s = \ 9ij C^, p# + J2 m - <% r - • z a + So, 



2 

where Sq is a function dependent on (q, z) only. Then, this function is not involved by the 
Gibbs-Appell equations and S can be replaced by 

S* = \ 9ij + ^ m v dijr v ■ d k v v z a . (6.1) 

V 

Now we show that this new function S 1 * assumes a very interesting expression. Let us introduce 
the functions 

djk(q) = Y m u dijv v • d k r u . 

Since 

= ^Z m v di{djr u • d k r u ) -^rn v (djr u • d ik r„) = dig jk - £ ikj , 

V V 

by a cyclic permutation of the indices we get 

Cijk £,ikj = 9i9jki (,jki ~\~ £,jik = 9jg k i, £,kij £,kji = 9 k gij. 

By summing the first two equations and subtracting the third one, since £ij k is symmetric in 
the first two indices, we get ^ k = Tij k , where 

Fij k = - (digj k + djg ki - d k gij) 
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are the Christoffel symbols of the metric tensor g^ (the coefficients of the Levi-Civita connec- 
tion). As a consequence, if we recall the definition (|4,3|) of G a p, the function 5* (|6.1|) can be 
written as 

S* = \ G a p z a z? + T ijk z a , 

and the Gibbs-Appell equations (|5.2|) assume the form 

G aP zP + T ijk ^^ k a = A a . (6.2) 
Then we can prove 

Theorem 6.1. The Gibbs-Appell equations t\5.2\i are equivalent to equations 

i a = G^(^-r jjfe W5), (6.3) 
where [G a/3 ] the inverse matrix of [G a p\ . 

Proof. Indeed, as we remarked in Section 01 the matrix [G a p\ is regular. If we apply the inverse 
matrix [G a ^] to equations (|6.2p . then we get equations (|6.3p , ■ 

Equations (|6.3p are the explicit form (or normal form) of the Gibbs-Appell equations (|5.2p . 



7 The dynamical equations of the first kind 



By setting 



z a 



dz a 
~dT' 



equations (|6.3p together with the constraint equations (|2.ip build up a first-order differential 
system, in normal form, in the unknown functions q l (t) and z a (t): 



d<f 

~dt 

dz a 

~dt 



Hence, we can summarize the results so far obtained in the following 



(7.1) 



Theorem 7.1. Let Q be the configuration n-manifold of a mechanical system, with local La- 
grangian coordinates q = (q l ). Let TQ be the tangent bundle of Q with canonical coordinates 
(q,q) = (q l ,q 1 )- Let C C TQ be a submanifold representing kinematical time-independent con- 
straints. Assume that the constraint submanifold C is locally described by parametric equations 
q l = ip l (q,z) with m < n parameters z = (z a ), such that the n x m matrix 



dtp 1 
~dz^ 



has maximal rank m. If the constraints are ideal, then the actual motions are represented 
by functions q l = q l (t) satisfying the differential system (7.1 1, where are the Christoffel 



symbols of the metric tensor gij associated with the kinetic energy K = ^ gij q 1 ^ o,nd A a = A% ^ 
where Ai are the Lagrangian active forces. 



a ■ 



3 Any solution of the differential system (|7.ip is of course a set of functions q l {t) and z a (t). But, after the 
integration, we can get rid of the functions z a (t). 
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Remark 7.1. The dynamical system (|7.ip is the first-order system associated with the vector 
field Z on the constraint submanifold C, whose components (Z l ,Z a ) w.r.to the coordinates 
(q, z) are given by 

Z a = G^{A p -T ijk ^^^). (7.2) 

The Z a are the 'vertical components' of Z w.r.to the projection onto Q. The actual motions are 
the projections onto Q of the integral curves of Z. 

Remark 7.2. This theorem provides a first ' recipe 1 for writing the dynamical equations for 
non-holonomic systems with linear or non-linear ideal constraints: 

1. Choose Lagrangian coordinates (q 1 ), write the kinetic energy of the system K = ^ Qij q l q^ , 
and extract the n x n matrix [gij]. 

2. Choose parametric equations q 1 = ip l (q,z) of the constraints, compute the m x n ma- 
trix [ipa], and check its rank. If the constraints are initially expressed by implicit equations, 
then use (for instance) the method illustrated in Remark 17.31 below for finding parametric 
equations. 

3. Compute the m x m matrix [G a p] = [gijip l a i^a\, and the inverse matrix [G a/3 ] = [G Q/ 3] _1 EI 

4. Write the Lagrange equations of the free motions (i.e., with only active forces Ai) in the 
form 

gijij* = Li(q,q) (7.3) 
(note that the formal expression of Lj is L, t = Ai(q, q) — T^ki q h q k ) and compute 

Zi{q, z) = Li{q, = Mq, ^) - T hkl (q) ^V- (7-4) 

5. Compute Z a = Z^i and Z a = G a/3 Zp. 

6. Write the dynamical system 

^ = Z"(q,z). (7.5) 

Remark 7.3. When the constraint submanifold C C TQ is described by a system of implicit 
independent equations of the kind 

C a (q,q)=0, a = l,...,r, r = n — m, (7.6) 

then we have to transform these equations into parametric equations. The choice of the pa- 
rameters z a is completely free and it is only a matter of convenience, depending on the explicit 
concrete form of the dynamical equations (|7.5j) we get. Anyway, since the matrix [ip l a ] does not 
have the maximal rank, equations (|7.6p can be solved w.r.to m of the n Lagrangian velocities q l , 
say - up to a reordering - w.r.to q a , a = 1, . . . , m. This process leads to considering as param- 
eters these Lagrangian velocities: z a = q a . It works very well for linear or affine constraints, 
where equations (|7.6j) have the form 



Cf(q)q i -b i (q)=0. 

Another possible choice of the parameters, for the linear constraints only, is that related to the 
use of quasi-velocities or quasi- coordinates - see for instance [9]. 



4 Note that according to this recipe we do not have to compute the inverse n x n-matrix of [g%j], but only the 
inverse of [G a ^], whose dimension is m < n. 
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Remark 7.4. For the analysis of the qualitative (or quantitative) behavior of a non-holonomic 
mechanical system (like stability equilibrium states, small oscillations, numerical integration, 
etc.) we can apply to Z all the known theorems about dynamical systems. For instance, a first 
integral is a function F(q, z) such that 

^,z)^ + Z-(q,z)^ = 0. (7.7) 



Remark 7.5. The singular points of the dynamical system (|7.ip are the solutions (q, z) of the 
simultaneous equations 

Z a (q,z) = Ziii = Q. (7-8) 

In the case of a homogeneous quadratic constraint, ip z (q,z) = \ ^ l a p{q) z a z@, equations (|7.8j) 
become 

^(q)z a z? = 0, 

Zi^zP = 0, (7.9) 

being = ip l aj3 (q) z@. This shows that, whatever q and Zi, singular points are given by z a = 0. 
But for z a = the matrix [ift^] does not have the maximal rank, since all its elements vanish. 
Hence, at these singular points the constraint C is not regular. 

Remark 7.6. The geometrical picture of the above results gives an intrinsic meaning of the 
objects we have introduced. Any vector v £ TC can be represented by a sum 

v = v i d i + v a d a , 

where di = d/dq l and d a = d/dz a are interpreted as pointwise independent vector fields on TC: 
at each point x of C they span the tangent space T X C. The vectors d a are 'vertical' i.e., they are 
tangent to the fibers F q of C . Hence, we call v a the vertical components of v, while we call v l 
the basic components. For instance, the basic components of Z are Z % = ip l (q,z) and the 
vertical components are Z a . The functions ip^ have the role of transforming Latin components, 

labeled by Latin indices h, k, . . ., into Greek components, labeled by Greek indices a, (3, 

For instance, when it is applied to a one- form (covariant vector) Zidq 1 , we get a vertical one 
form Z a dz a , and when it is applied to the covariant metric tensor gij , we get a metric tensor G a p 
on the fibers of C, so that, by raising the indices of Z a by the contravariant metric G al3 , we get 
a vertical vector field Z a d a . 

Note 7.1. I did not find equations (17. ip in the recent and old articles I have consulted. In fact, 
it is rather surprising that the simple idea of considering the parametric representation of the 
constraints as a part of the dynamical equations does not appear in the major textbooks and 
treatises on non-holonomic mechanics. Only recently this idea appeared in a paper of Massa and 
Pagani [8]. Their general approach, which is based on the jet-bundle theory and deals with time- 
dependent constraints, leads to the introduction of the vector field Z. The elementary approach 
presented here is of course quite different and leads, for instance, to different expressions of the 
vertical part of Z. Our equations (|7.ip should be compared with equations (3.5b) and (3.15) 
of [8]. The difference is that the second equation (17. 1|) is written in terms of the Euclidean 
vectors F^, v u , while equation (3.15) of [8] Z is written in terms of the Lagrange equations, but 
still in an implicit form. 
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Remark 7.7. For linear constraints, the approach presented here is more general than that 
of Caplygin - see [9], Ch. Ill, § 3, where the coordinates q l are divided into two groups, say 
(q a , q a ), with a = 1, . . . , m and a = m + 1, . . . ,n. The constraint equations are assumed to be 
of the form 

a 

where the coefficients and the Lagrangian L are assumed to be independent from the coor- 
dinates (g a )H 



8 The dynamical equations of the second kind 

About the method for writing the dynamical equations of a non-holonomic mechanical system 
so far illustrated two remarks are in order: 

• It lies on a parametric representation of the constraint C (however, the vector field Z does 
not depend on the chosen parametrization) . 

• It does not give any information about the reactive forces. 

Here, we propose an alternative method for writing the dynamical equations which is based 
on any implicit representation of C by a system of independent equations 

C a (q,q)=0, 

and which provides a way for evaluating the reactive forces. 

The Lagrange equations (|3. 14[) - see Remark 13.41 - are equivalent to the dynamical system 




=h= -T l hk q h q k +A l + I? x , 



on the tangent bundle TQ of the configuration manifold Q. Here, A 1 = g %3 Aj and R l x = g lJ X a Cj 
are the contravariant components of the Lagrangian active and reactive forces, respectively. 
The label A points out that the reactive Lagrangian forces depend on the a priori unknown 
Lagrangian multipliers A = (A a ), according to Theorem 13.11 and equation (|3.11[) . 

For a better understanding of what we are going to do, it is useful to consider the following 

Definition 8.1. We say that a vector V on TQ is vertical if it is tangent, at each point where 
it is defined, to the corresponding fiber of TQ. This is equivalent to say that it has the form 

V-V£. (8.2) 



As a consequence, the vertical part of the vector is given by 



d_ 

dq l 



5 In the history of the non-holonomic systems we can find the famous equations of Maggi, Volterra, Voronec 
and Caplygin, dealing with linear constraints. The comparison of these equations with our approach is left to the 
reader, who can find a detailed discussion in in [5], Ch. 3. A neat illustration of the non-holonomic dynamical 
equations, with the essential classical and recent bibliography, can be found in the book [2]. 
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where the active and reactive forces are represented by the vertical vectors 

A = A 1 — — r, R\ = R\ — — 7. 

dq 1 A dq 1 ' 

Hence, the vector X A is decomposed into the sum 

X A = X G + A + R A , (8.3) 
where X G represents the geodesic flow, 

f — = <f 

X G = * (8.4) 
( dt ~ ~ Thk<1 Q ■ 

We can consider X A = X G + A + R A as a family of vector fields, depending on the Lagrangian 
multipliers. However, it is a remarkable fact that we can obtain an explicit form of them, as 
functions of the kinematical states (q, q) only. 

Theorem 8.1. Let [G ab ] be the symmetric matrix defined by 

G ab = g ij C?C b . (8.5) 

Let [G a b] = [G *]" 1 be its inverse, and C a% = Cj. If the constraints are ideal, then the 
Lagrangian multipliers and the Lagrangian reactive forces are well determined functions of(q,q): 

X a (q, q) = G ab (Cf q h q k ~ A*) - cf d^) , (8.6) 
R\q, q) = G ab C ai (c* (T{ k q h q k - A*) - cf d,C b ) . (8.7) 

Proof. First of all, we observe that the matrix [G ab ] is regular, since the vector fields C a are 
independent. Hence, the inverse matrix [Gab] is well defined. In order to satisfy the constraints, 
the vector field X A must be tangent to the constraint submanifold C. This condition is expressed 
by equations 

(X A ,dC a )=0, (8.8) 
to be satisfied at least on C. In components, these equations read 
q i diC - + (A 1 - V hk q h q k + &) Cf = 0, 

i.e., 

R'Cf = Ct(Ti k q h q k -A i )-q i d i C a . 

Note that the right hand side 

A a (9, q) = C? (r hk q h q k - A 1 ) - cf diC a (8.9) 

is a known function of (q, q). Then equation R l Cf = A a assumes the form X b Cf C b cf J = A a , i.e., 
A;, G ab = A a . By applying the inverse matrix [G a b] we get equation (|8.6p and equation (|8.7p . ■ 



The explicit form (|8.7p of the reactive forces allows us to state 
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Theorem 8.2. The actual motions of a mechanical system with regular and ideal non-holonomic 
constraints represented by a submanifold C C TQ are the integral curves based on C of the vector 
field 

D = X G + A + R, (8.10) 

where the components R 1 of the vertical vector R are defined by t \8. 7[ ) . If we introduce the symbols 

^ = G ah C ai C b \ ("" ■ Ci = G ab C b \ (8.11) 

then the explicit expressions of the first-order differential system associated with D and of the 
reactive forces are 



D 



dq l 

~dt 

da 1 



1.12) 



and 



# = 7T ij ' (T hkj q h q k - Aj) - qi d,C a Cj, 
respectively. 



.13) 



Remark 8.1. This last theorem provides a second 'recipe' for writing the dynamical equations 
for non-holonomic systems with linear or non-linear ideal constraints: 

1. Choose Lagrangian coordinates (q l ), write the kinetic energy of the system K = ^ gij q l tf , 
extract the n x n matrix [gij], and compute the inverse matrix [g l3 \. 

2. Take the constraint equations C a (q,q) = and compute, in the order, the following ma- 
trices: 



[Cf] 

[C ai \ 
[G ab 

[Gab 
[CI] 



dC a 



(seek the singular states), 



7T 



ye" ■ 

= [G ba ] = g * (ye 1 ; = [C ai C b ], 
-- [G ab ]~ 1 , 
[G ab C% 



7T 



J'] 



[c ai ci], 



7T ' 



dC a 
dq % 



[g- 

[d t c a ] = 

[djC a Ci\. 

3. Write the Lagrange equations for the free motions in the form 
gij q j = Ai-T hki q h q k , 
and keep in evidence the functions 
Li(q,q) = Ai-T hki q h q k . 
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4. Compute the vector 

I)' = [<,'■> - TT lj ) Lj - q> i)jC" (•'„. 

5. Write the differential system (|8.12p . 



D 



~dt 
dcf 
V ~dt 



D\ 



6. Its solutions q l (t), q l {t), with initial conditions belonging to C, describe the actual motions 
of the system. 

7. Evaluate the reactive forces along any actual motion by means of equation (|8.13p . 



n ij Lj - q j djC a C l a . 



$.14) 



Remark 8.2. The case of single constraint equation C(q,q) = 0. In this case the above- 
given recipe can be applied by setting a = b = 1. Items 1 and 2 of the general recipe still hold. 
However, since some of the above matrices reduces to scalar functions or to vectors, the index 1 
can be omitted or replaced by *: 

[&] ± \g*t Cj], 
d = G' 1 C\ 

yo\ = [c i ci] = G~ 1 [c*c J 'j, 



[9 J -7T J J, 

[djCCt]. 



Then, follows items 3-7 of the general recipe. 



9 Illustrative examples 

As shown above, for writing the dynamical equations of a non-holonomic system we can apply 
two methods: the first method is established by Theorem 17. II and the corresponding recipe is il- 
lustrated in Remark l7.2t the second method is established by Theorem l8.2l and the corresponding 
recipe is illustrated in Remark 18.11 

Let us see how these two methods work by concrete examples. We begin with two paradig- 
matic and simple examples of linear non-holonomic constraints, the 1 skate'' and the 'vertical 
rolling disc'. Then, we shall consider two more demanding examples: l two co-axial rolling discs' 
and l two points with parallel velocities' . This last one is a genuine non-linear non-holonomic 
system. 

In illustrating examples of application of a theory it is not customary, in general, to provide 
detailed calculations - which usually are left to the reader. Here, however, it is worthwhile to 
disregard such a custom in order to compare the effectiveness of the two methods (mainly the 
length of the calculations) applied to a same mechanical system. 
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9.1 The skate 



This mechanical system is made of a homogeneous rod (material segment) sliding without friction 
on a plana^t The configurations of the skate are determined by the Cartesian coordinates (x, y) 
of the center of mass (i.e., of the segment) G and by the angle 6 of the rod w.r.to the x- 
axis. The configuration manifold Q is R x Si and natural ordered Lagrangian coordinates are 
{q 1 -, Q 2 -, Q 3 ) = ( X ,U,0)- The velocity vg = [x,y] of the mass-center is constrained to be parallel 
to the rod. This constraint is then represented by a single equation: 



x sin 6 — y cos 6 = 0. 



(9.1) 



The kinetic energy is given by K = i m (x 2 + y 2 ) + \ l 6, where m and / are the mass and the 
moment of inertia w.r.to G (i.e., w.r.to the line orthogonal to the plane through G), respectively. 

(i) First method. Since dim(TQ) = 6 and dim(C) = 3, we need two parameters (z l ,z 2 ) for 
describing C. We can consider the parametric equations 



x = z cost/, y = z suit/, t> = z . 
Then, we compute the necessary matrices and vectors: 



\90\ 



m 



\z v cos l 




m 
I 

i 



Z Sill! 



COS I 





8111 ( 







1 



This matrix has maximal rank everywhere: the constraint is regular. 



[Gap] 



m 
/ 



\G 



aP] 



t 

7 



Since, 



dK 



\mx 



my , 16], 



dK 

dq l 



the Lagrange equations for the free motions g\j q 3 

m'x = A\, my = A2, 1 6 = A3. 
Hence, L« = Ai(q,q), Z { = Ai(q,ip) and 



[0,0, 0], 

Ai - T hki q h q k read 



[Z a ] = [i> a Z i ] 



A\ cos 6 + A2 sin ( 

^3 



The dynamical equations are 



dx 

~dt 
dy 

dt 
<W_ 

~dl 



Z COS I 

z l sin t 



dz 1 / A n 

= — (Ai cos 6 + A 2 sin I 

dt m 

dz 2 



dt 



Az 
I ' 



e Quite similar classical examples are that of two material points linked by a massless rigid segment [5], p. 23 
and 63, and the Caplygin sleigh [S], Ch. Ill, § 3, Examples 2 & 5, and Ch. V, § 4. Another example of this kind 
is examined in [4], § 4.1 & § 4.2. 
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(ii) Second method, for a single constraint equation - Remark 
^-0 



C = x sin 8 — y cos 8 = 0. 



[9 



[Ci} = 






dC 



,a 
0} 



dq 

\C l ] = [g ij Cj] 



[sin 8 , — cos 8 , 0] 



sm( 



rn 



cos ( 



, 



G — C l Ci — — , 
m 

[CI] = G-^C*] = [sin^ , - cos 8 , 0] = [Q], 

sin 2 8 — sin cos 8 
— sin 9 cos cos 2 


cos 2 8 sin 8 cos 
sin 8 cos sin 2 8 




[tt^'I = G~ l \C l C j ] = — 
m 



IfJ 



m 



m 
i 



[diC] 



dC 



dq 1 



[0 , , x cos 8 + y sin 0] , 



[q j djC C*] = 8 (± cos 8 + y sin 0) [Cj] . 



Since 



dK 



dq 1 



[m x , my , 1 9], 



dK 

dq 1 



[0,0, 0], 



the Lagrange equations read 
mx = Ai, my = A2, 



18 = A,. 



Hence, Lj = A{. We have all the ingredients for computing the vector D l 



D 



D 



cos 8 

m 
sinf 



rn 



{A\ cos 9 + A2 sin 8) — 8 (x cos 8 + y sin #) sin 6>, 
{A\ cos 9 + A2 sin 0) + (± cos # + y sin 0) cos 9, 



D 3 = jA 3 



and the differential system (|8,12p reads 



dx 
dy 

-dt =V ' 
d8 _ . 



(ii; cos ( 



<iy sin 6 



(it m 
dt 6 



(A\ cos 8 + A2 sin 9) — 9 (x cos # + y sin #) sin 0, 
(j4i cos 6 + A2 sin 0) + 8 (x cos + y sin 0) cos 6*, 
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9.2 The vertical rolling disc 

A material disc of radius R running on a plane is kept perpendicular to it by massless and 
frictionless devices. The configuration manifold is Q4 = M 2 x Si x Si, with coordinates 
(g 1 , g 2 , g 3 , g ) = (x,y,9,ip), where (x,y) are Cartesian coordinates of the center P of the disc 
(i.e., of the point C in contact with the plane), 9 a rotation angle of the disk around its axis, 
and ip an angle giving the orientation of the axis (see Figure 4, with 8 = 81). Constraint: the 
disc rolls on the plane without sliding. Let (i,j,k) be the unitary vectors associated with the 
(x, y, z)-axes. The unitary vector k is associated with the oriented angle ip. Let u be the unitary 
vector associated with the oriented angle 9. Then, u = cos ip j — sin^ i. The angular velocity 
u> is given by u) = 9 u + if) k. The velocity vc of the point C is given by vc = vp + u x PC, 
where PC = — R k. Hence, 



vc = ii + yj-(^u + ^k) xi?k = xi + yj 
Since u x k = cos ip j x k — sin ip i x k, we get 



R 9u x k. 



vc = xi + yj — R9 (cos ip j x k — sin -0ixk) = xi + yj — R9 (cos ip i + sin ip j) 
= (x — R9 cos + (y — R9 sin tjj) j. 

The kinematical constraint = is then represented by the following two linear equations 



C 1 = x 
C 2 = y 



R cos ip 9 = 0, 
R simp 9 = 0. 



(9.2) 



(i) First method. Assume that the center of mass of the disc coincides with its geometrical 
center. Then the kinetic energy is given by 



K = \ mix 2 + y 2 ) + \ (A9 2 + B1J; 2 ), 



(9.3) 



where m is the mass, A and B are the moments of inertia w.r.to the axis of rotation and 
a diameter, respectively. Thus, 



[9ij] 



m 

m 

A 

B 



From the constraint equations (|9.2I) we get the parametric equations 



x = R cos if) z 1 , 
y = R simp z , 



if) 



(9.4) 



Thus, 



w 

[4 



[R cos ip z 1 , R sin ip z 1 , 

i? cos if) R sin -0 10 
1 



(a = 1,2, index of line). 



This matrix has maximal rank, thus the constraint is regular. It follows that 

G aP = giji&ijPp = m ip x J>\ + m iP 2 JPp + A + B 

Gn = mip\ip\ + mipltp 2 + Atpfipf + B ipfipf = mR 2 cos 2 if) + mR 2 sin 2 ip + A = mR 2 + A, 
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C22 - 
fG Q/3 l 



mi? 2 + A 



0, 

B, 







5 



1 



mi? 2 + j4 





1 

B 



Since, 



dK 



dq l 



[mx , my , j4 , Bi/j], 



dK 
dq l 



[0,0,0,0], 



the Lagrange equations for the free motions gij <p = Ai — q h q k read 

mx = A\, my = A 2 , A9 = A3, Bip = A 4 . 
This shows that Li = A{. Hence, Zi = Ai(q,tjj), and 



\Z, 



[Z a ] = [G Q/3 Z p \ 



At R cos tp + A 2 R sin ip + A 3 
A 4 

A x R cos ip + A 2 R sin i> + A 3 
mR 2 + A 
JU 
B 



Thus, the dynamical equations are 
dx 



dt 
dy 

dt 
d9 

lit 

d4 



R cos ip z , 
R sin ^2; , 



dz 1 
~dt 
d£_ 
dt 



(mR 2 + A)' 1 (Ax R cos tp + A 2 R simp + A 3 ), 
B- X A 4 . 



(9.5) 



(ii) Second method. Recall the constraint equations (|9.2p . Then, 
[Cf] 

[C ai 



~dC a ~ 




10 — R cos if; 


" 


dq l 




1 — R sin ^ 


_ 



ab] 



\/ j cf] 



C ai c b 



i 

m 

h 



4 cosV' 



sin^ 



m + "X cos2 ^ "X sm ^ cos ^ 
A- sin ^ cos ^7 + "X sm2 ^ 



G = det[G ab ] 

[G a b] = G 1 



1 



+ 



A 

R 2 



G- 1 



A 

m 2 A 



m* mA ' mR? + A 

— + ^ sin 2 ^ — "X s i n V 7 cos * 

— V sm ip cos V — + -j- cos 
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[C l a ] = [G ab C bl ]=G 



±(± + ^sin 2 iP) -^sin^cosV *Cf 
-^sin^cos^ X(i. + ^ cos 2^) +C 3 o 



where 



Hence, 



1 i? 2 



m ^4 
i? 



R 



*C 1 = I 1 — — sin ip I — — cos-i/ 1 + I t- sin ^ cos ^ I — - sin?/> 



,4 



R 2 



A 



R 



A 



R 3 



R 3 



mA 
A 

A 



cos ip ttt sin ip cos ^ + —rz sin ?/> cos ib 



A 2 

sin ^ cos ip 
sin ?/> cos ^ 



A 2 



R 

mA 



cos 0, 



R 
A 
R 
A 



1 


+ 


R 2 


m 


A 


1 


+ 


R 2 


m 


A 



R 
A 
R 
A 



sin ip 
sin ip 



R 

mA 



simp. 



[Cj] = G^ 1 



^(^ + l sin2 V') -^sin^cosV -^cosV 
-^sinVcos^ £ + f cos 2 V) "^sinV 



Let us compute vr*-? = C a * Ca = vr Jt . Let us set *ir t: > = Gif'K Then, 



[ir ij ] = CT 1 



£j sin ^ cos ^ *7r 14 



m 2 A 



-Sfj sin ib cos t/? A- (— + 4r cos 2 -0) *7r 24 



*7T 



31 



*7T 



41 



*7T 



*7T 



,1 
32 

42 



*7T 33 





where 



*7T 



13 



*7T 
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1/1 A • 2 ,' 

— I h — sin ip 

m \m A 

i/i R 2 • 2 , 

— I 1- — sin ip 

m \m A 



RP_ 
mA 
R 
m 2 A 



sin ib cos ip 



sin^. 



It follows that 



~ cosip 



R 

-— cosip ) H 



ii 2 . . 

sin cost, 

mA 

sin 2 ip cos f/j 



A 



cosip + 



mA 2 
m \m A 



cos 2 ^ 



i? 


sinib 


A 


R 


cos ip 


m 2 A 


' R 


sin ip 


A 



m z \ m A r ' 



sin ip cosip 



^sinVcosV ^Q + TT cos 2 V) 



TJ COS^ 



m 2 A 





— tt sin^ 

m z A ™ 



-jl COS^ 
m^A r 



fj sin ip 





R 2 
mA 2 
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[ g V - rfi] = G- 1 



— *t-t sin ib cos ib 



sin ^ cos -0 



R 2 



R 2 ./. „™ ./. G _ 1 / J_ , ros 2 /A _R 



5-^ suit/; 




cos ?/; 



sin-0 



fl 2 



G 



771 



mi? 2 + A 



G 



;£ 2-r sin t/j cos ^ 



m m' 3 v m 
?2 



m' 



sin ^ cos ^ 



^4 cost/; 



L (£ + f cosV) ^sin^ 



m 2 A 
G R 2 



1 



R 2 + A 



111 





ii 2 cos 2 ib R 2 sin ^ cos ip R cos ^ 

R 2 sin ?/> cos ^ R 2 sin 2 ^ sin ip 



A mA 2 





G 

B 



R cos ^ i? sin ip 


Recall once more equations (|9.2p . Then, 



1 









iR 2 +A 
B 



[diC° 



dC a 
dq l 



R simp 9 
- R cos ip 9 



[q i d l C a 



R smip9ip 
- R cos ipOip 



Let us set X i = q j djC a C*. Then, 

1/1 R 2 



GX 1 = (R sinipdip) 



+ ^- sin 2 ^ ] ] + (--R cos ip9ip) 
m \ m A 



mA 



sin ip cos 7/^ 



R . 1 „ ; R 1 a ; R 

— - sm tp9ip -\ sm ip 9 ip = — = 

m z mA m l 

R (mR 2 + A) 



sin ip 9 ip ( 1 + 



mR 2 
~A~ 



in 



2 A 



simp 8 ip, 



GX 2 = (R sinip9ip) 



mA 



sin ip cos ip + (-R cosip9ip) — 1 — — cos ip 

m V 777 A 



1/1 i? 2 



2 COS IpUtp - COS Iptitp ~ 



77? 



GX 3 = (i? sinipdip) 



mA 
R 

mA 



m 2 A 

cos ip ) + (-R cos 7/; d ip) 



cos "0^7/^ 
i? 



sin 1 !/; = 0, 



GX 4 = (R simp9ip)(0) + (-12 cos^0^)(O) = 0. 



i? {mR 2 + A) 
to 2 ^4 



sin ip9ip 



R (mR 2 + A) 
m 2 A 



cos ip9 tp , , 



^ [sin ^ , — cos tp , , 0] 
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Now we are able to compute the components D l = [g l] — ir 13 ) Lj — q J djC a C % a of the vector D: 

D 1 = — -3r r (Ai R 2 cos 2 ip + A 2 R 2 sin ip cos ib + A 3 R cos ip) - R 6 ip sin ip 

mR z + A 

Fi, cos ib 

(A 1 R cos ip + A 2 R sin tb + A 3 ) - R 9 ip sin ip, 



mR 2 + A 
1 



D 2 = — -= (Ai R 2 sin ib cos ib + A 2 R 2 sin 2 ib + A 3 R sin ib) + R 9 ib cos ib 

mR z + A y ' 

R 'smib 



m 



R 2 + A 



(Ai R cos ib + A 2 R sin ib + A 3 ) + R 9 ib cos ip, 



, 1 

-D = —-5 T iii? cosV + A 2 .R sinV + ^3) 

mR z + ^4 

(A x i? cos ib + A 2 -R sin ib + A 3 ) , 



mR? + A 

- *± 
B ' 

The resulting dynamical system (|8.12p is 



D 



dx 




dx 


R cos ip 


~dt 


= x, 


~dt ~ 


mR 2 + A 


dy 




dy _ 


R sin-i/j 


dt 


= y, 


dt ~ 


mR 2 + A 


dd 


= 9\ 


d9 


1 


~dl 




dt ~ 


mR 2 + A 


dip 


= ip\ 


dip _ 


A 4 


~dt 


~dt ~ 


B ' 



(Ai R cos ib + A 2 R sin ip + A 3 ) - R 9 ip sin ip, 
(Ai R cos ip + A 2 R sin ip + A 3 ) + R 9ip cos ip, 

(A 1 R cos ip + A 2 R sin ip + A 3 ) , 



By introducing the new variables 

__ . mR 2 + A ^ . mR 2 + A . 2 

X = x, Y = y, ® = {mR 2 + A)9, 

it assumes the more compact form 

— = x -r- = cos ib (A x R cos ip + A 2 R sin ip + A 3 ) - 9 ip sin ip, 
dt ' dt 

clY dY 

— — = Y, — — = sin ip (Ai R cos ip + A 2 R sin ip + A3) + ip cos V>, 
at at 

dQ ■ dQ 

~df = ®> — = (AiR cosip + A 2 R sin^ + A3) , 

d f = i,, d± = M 

dt dt B' 

Note that in these new variables, by considering also 

Z 1 = (mR 2 + A)z l , 

the differential system (19.51) obtained by the first method reads 

dX ! 

— = cos ib Z , 

dY dZ 1 

—— = sin ip Z , —— = A\R cos ip + A 2 R sin ip + A3, 

at at 

dt dt 
d^ _ 2 
dt ' 



(9.6) 



(9.7) 



(9.8) 
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The two systems (|9.7p and (|9.8[) are in perfect agreement. 

The above detailed calculations show that for the rolling disc the first method is much shorter 
than the second one. 

9.3 Two co-axial rolling discs 

Two identical material discs of radius R running on a plain are joined by a massless common 
axis, along with they can slide without friction. The configuration manifold is Qe = ]R 2 x Si x 
Si x Si x K, with Lagrangian coordinates (q , q 2 , q 3 , q 4 , q 5 , q 6 ) = (x,y,6i,02,ip,a), where (x,y) 
are Cartesian coordinates of the center P\ of one of the two discs, 8± and 02 are the angles of 
rotations around the common axis, ip is the angle giving the orientation of the axis, and a is the 
distance between the centers (see Fig. 4). 




Figure 4. Co-axial rolling discs. 

Constraint: the discs roll on the plane without sliding. For each disc this constraint is 
represented by linear equations of the kind (|9.2|) . 

X\ — R COS Ip Q\ = 0, X2 — R cos 1p02 = 0, 

(9.9) 

yi — R sin ip 0\ = 0, V2 — R sin ip 62 = 0. 

However, the coordinates of the two centers are related by equations 

x 2 = xi + a simp, y 2 = y\ — a cos ip. 
By differentiating these equations we get the link between the velocities, 

±2 = x\ + a cos ip ip + sintp a, ij2 = yi + a sinip tp — cos ip a. (9.10) 

By inserting these relations into equations (|9.9p . with x\ = x and y± = y, we get the final 
constraint equations 

x — R cosi/jOi = 0, 
y — R s'mtb 9\ = 0, 

. (9.11) 

x + a cos ip ip + sin ip a — R cos ^#2 = 0, 

y + a sin tp ip — cos ip a — R sin ip 62 = 0. 

Since for a single rolling disc the first method is faster, we limit ourselves to apply the first 
method to the case of two discs. 
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Equations (|9.1ip show that the constraint submanifold C C TQq has dimension 12 — 4 = 8. 
Hence, we need two parameters z a for the parametric representation. Let us choose z 1 = 6\ and 
z 2 = #2- This means to solve the linear system (|9.1ip w.r.to (x,y,ip,a). The result is 



X = 


R cos ip z , 


d = 0, 


y = 


R sinip z , 


9i = z 1 


i> = 


— (z - z ), 


9 2 = z 2 



Equation a = 0, a = constant, exhibits the intuitive fact that, under the pure-rolling con- 
dition, the distance a between the two discs remains constant. Hence, we can reduce the 
configuration manifold Qq to Q5 = M 2 x Si x Si x Si, with coordinates (q 1 , q 2 , q 3 , q 4 , q 5 ) = 



(9.12) 





0i,t 


h,ip)- The last equations reduce to 




X - 


= R 


COS')/' z 1 , 

sin if) z 1 , 8i = z , 




y = 


-- R 






R 

a 




with 


a = 


constant. Then, 




['J' 1 


= 


R 

R cos ip z 1 , R sinip z , — (z 2 — 

a 



and 



dip 




R cos ip R sin ip 


dz a _ 








R 

a 





1 



(9.13) 



This matrix has maximal rank, thus the constraint submanifold C is regular. The kinetic energy 
of the system is the sum of the kinetic energies of the two discs. According to equations ()9.3[) 
and (l9~T0jl . 

K=\ m(x\ + y\) + \ (AO 2 + Bij; 2 ) + \ m(x\ + y%) + \ {A9 2 2 + Bijj 2 ) 
= \ m(x 2 + yj + ± 2 2 + y\) + \ (AO 2 + Bij; 2 ) + \ {Ad 2 + Bij; 2 ) 
= \m( K x 2 1 + yl + {xi+a cos ip ip) 2 + (yi + a simp ip) 2 ^j + | (<9 2 + 0\) + Bip 2 

= h m ( x\ + yf + x\ + a 2 cos 2 ip ip 2 + 2a x\ cos ip ip + ij\ + a 2 sin 2 ip ip 2 



Thus, 



+ 2a Vl smipip) + A (Ql + Ql) + B tf 



\ m \2x\ + 2y\ + a 2 ip 2 + 2a ip {x\ cos ip + y\ sin ip)) + ^ (9f + 9%) + Bip 



m (x 2 + y 2 ) + ^ 01 + 9%) + (±ma 2 + B)ip 2 + maip (x cosip + y sin ip). 



[9 



2m 














2m 














.4 














.4 


cos ip 


22? sin^ 









^ cos^ 
^ sinV' 
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Moreover, 



G a /3 = 9ij C % = 911 i>\ i>\ + 522 i> 2 a i>l + 933 V$ + 544 V'q + 555 

+ 2g 12 ^ i>\ + 2g 13 ^ ^ + 2ffM ^4 + 2ffi5 ^1 ^5 + 2g ^ ^2 ^3 + 2g ^ ^2 ^4 
+ 2 525 i) 2 a $\ + 2 5 34 ^ + 2 5 35 ^ ^ + 2 5 45 4>t $\ 

= 2m V„ ^ + 2m^ ^| + $ + A ^ + (™a 2 + 25) ^ ^ 

+ ma COS ijj V'fl + ma ^ V>a V'fl; 

Gn = 2m ^ + 2m ^i ^i + A ipf V? + A ipj ^ + (ma 2 + 2J3) ^ ^i 
+ ma cos tpf + ma sin ?/> -^i V'i 



22 



o2 / R 2 \ 

2m i? 2 cos 2 ^ + 2m R 2 sin 2 ^ + A + A = 2mR 2 + A ( 1 -\ T ), 

a 2 \ a 2 J 

2m ip\ if>l + 2m iplipl + A ip% t/% + A tfc ip 4 + ( ma2 + 2B ) 



+ ma cos ip tp2 + ma si n ^ ^2 V'f = A — — + ma z + 25, 



i? 2 



G12 = 2m ipl ifi +2mi/}'i^ + Atpl^ + A ip\ + (ma 2 + 2B) ip\ ij% 



1.2 „/.2 



3 „/,3 



/,5 „/,5 



+ ma cos ip ip2 + ma s i n V' "01 V* 

R 2 R 2 
— A—^- + ma cos ^ i? cos ip + ma sin^ i? sin-0 = maR — A—^, 



2 „/,5 



and we obtain 
It follows that 



2mi? 2 + A (1 + 
mai? - AK 



maR - AK 



A K + ma 2 + 2B 



G = det[G a p] 



R 2 

2mR l + A ( 1 + 



A + ma 2 + 25 



maR — A—^r 

a z 



R 4 R 2 R 4 R 2 

2mA — + 2m 2 R 2 a 2 + AmR 2 B + A 2 — + ma 2 A + 2AB + A 2 — + ma 2 A— 



+ 2AB 



R 1 



2 2 D 2 

maR 



A 



2 R 4 



+ 2m^ 



R 3 



R 2 



[2mAR 2 + A 2 + ma 2 A + 2AB + 2m^ai? - m 2 a 4 + 2m 2 a 4 + 4ma 2 B) 



+ A (ma 2 + 2B) 
R 2 

= — (2mAR 2 + A 2 + ma 2 A + 2AB + 2m^ai? + m 2 a 4 + Ama 2 B) + A (ma 2 + 2B) 
R 2 

= — (m (2AR 2 + a 2 A + 2AaR + Aa 2 B + ma 4 ) + A 2 + 2AB) + A (ma 2 + 2B). 
a z 

We observe that the determinant G is a constant. The inverse matrix is 
A ^ + ma 2 + 2B A^r - maR 

\G afi \ = G" 1 

AZ? - maR 2mR 2 + A ( 1 + 
From the expression of the kinetic energy, 

K = m (x 2 + y 2 ) + ^ (^? + ^|) + (| ma 2 + 5) ■0 2 + ma^ (x cosip + y simp), 
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we obtain 

2mx + ma cos ip ip 
2my + ma sin ip ip 

A9 2 

(ma 2 + 2B)ip + ma(x cos ip + y sin ip) 








maip (y cos ip — x sin -(/>) 

The Lagrange equations for the free motions are 

2mx + ma cos ipip — ma sin ipip 2 = A\, 
2my + ma sin ip ip + ma cos ip ip 2 = A 2 , 
A(9i = A 3 , 
A0 2 = A*, 

(ma 2 + 2B)ip + ma(x cos ip + y sin ■;/>) — ma(x sin ip — y cos f/")^ 
= maip(y cos V> — isin^) + A5. 

They show that 

Li = + masimp ip 2 , 
L 2 = A 2 — ma cos ipip 2 , 
L 3 = A 3 , 
L4 = A4, 

L5 = A§ + mai/i (y cosV' — i sin^) + ma(x simp — y cos tp = A§. 

Thus, due to the parametric equations (|9.12p . 

Z x = A l + ^f. ( z 2 - z 1 ) 2 sin^, Z 2 = A 2 - tzf- (z 2 - z 1 ) 2 cos V, 
Z 3 = A 3 , Z 4 = At, Z b = A 5 . 

Let us compute Z a = Ziip l a - recall (|9.13p : 

Z x = (A 1 + ^f- (z 2 - z 1 ) 2 simp) R cos ip + (A 2 - ^f- {z 2 - z 1 ) 2 cos ip)R sin ip 
"#^3 + 4. 

= Ax R costp + A 2 R simp - ^ A 3 + A A = R(A X cosip + A 2 simp - ± A 3 ) + A 4 , 
Z 2 = %A 3 + A 5 . 



dK 

dq l 



'dK 
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Thus, the dynamical equations (|7.5p associated with the vector field Z are 
R cos ip z , 



R simp z , 
— {z - z L ), 



dz 1 

— = R(Ai cosip + A 2 sini/j - ± A 3 ) + A 4 , 
at 



z\ 



dz^ 
dt 



%A 3 + A 5 . 



dx 

~dt ' 

dx 

~dt ' 

d4 

~dt 
dBi 
~dt 

d02 _ 2 

dt ~ z ' 

9.4 Two points with parallel velocities 

Two material points Pi = (x±, yi) and P2 = (X2, 2/2) running on the Cartesian plane R 2 = (x, y) 
are constrained to have parallel vector- velocities vi and v 2 - This is an example of non- linear 
non-holonomic constraint, since it is expressed by the quadratic homogeneous equation 

C = x 1 y 2 -x 2 y 1 = 0. (9.14) 
The configuration manifold is Q4 = R 4 with ordered Lagrangian coordinates 

(<7\<? 2 ,<? 3 ,9 4 ) = (xi,yi,x 2 ,y 2 ). 
The kinetic energy is K = % mi {x\ + y\) + \ iri2 (i 2 + Hence, 



[9i, 



(i) First method. Since dim(Q) = 4 and dim(C) = 7, for a parametric representation of the 
constraint x\ y 2 —x 2 yi = we need three parameters (z , z 2 , z 3 ). Let us consider the parameters 
(z a ) = (z 1 , z 2 , z 3 ) = (p, cr, 0) and the parametric equations 



mi 














mi 














m 2 














m 2 



X\ = p COS 0, 

yi= p sin0, 



X2 = cr cos ( 
y 2 = a sin 



(9.15) 



The meaning of the parameters is the following: p = = x\ + yf, a = ^ 2 = %2 + i)2i an d Q is 
the angle of the two vector velocities w.r.to the a;- axis. Then we find: 



[Gap] = [9ij C^] 



COS V 8111 V 



- p sin 9 p cos 9 

mi 

m,2 




cos 9 sin 9 
a sin 9 a cos 1 






[G 



a/3] 



1 

— 

mi 

— 

m 2 





mip 2 + m 2 a 2 




1 



mi/O 2 + m2<T 2 



(q index of line), 



Dynamical Equations of Non-Holonomic Systems 



29 



Since, 



dK 



dq 1 



[mixi , miyi , m 2 x 2 , m 2 y 2 \, 



dK 
dq l 



[0,0,0,0], 



the Lagrange equations for the free motions gij g J = Ai — q h q k read 

m 1 x 1 = A 1 , m 1 y 1 = A 2 , m 2 x 2 = A 3 , m 2 y 2 = A A . 
They show that Zi = Ai. Hence, 



(9.16) 



[Za] = WM 



[Z a ] = [G al3 Z p ] 



Ax cos 6 + A 2 sin 9 
A3 cos 9 + A± sin 9 
p (A 2 cos 9 - Ax sin 9) + a {A A cos 9 - A 3 sin 9) 
A\ cos 9 + A 2 sin 9 
mi 

A3 cos 9 + A4 sin 9 
m 2 

p (A 2 cos 9 — Ax sin 9) + a (A± cos 9 — A3 sin I 
m\p 2 + m 2 a 2 



and the differential system associated with Z is 

^=pcos9, dp _ Ax cos 9 + A 2 sin9 



dyx 
dt 
dx 2 

~dt 

dm 

dt 



p smu, 
■ a cos 9, 
crsin#, 



dt mi 

da A3 cos 9 + A4 sin 9 

dt m 2 

d9 p (A 2 cos 9 — Ai sin 0) + a (A^ cos — A3 sin I 
eft ~ ~ 



(9.17) 



niip 2 + m 2 a 2 



where the Lagrangian active forces are in general known functions of (xi,yi,x 2 ,y 2 ) and (p, a, 9). 

In the special case of an inclined plane we have Ai = mi g, A3 = m 2 g, A 2 = A4 = 0, and 
the system (|9.17p becomes 



dxi 
~~dt 
dyi 
dt 
dx 2 

~dt 

dy 2 
dt 



■ p cos p, 
p sin0, 

: a cos 9, 
crsin#, 



dp . 
da 

-dI = 9C ° S( 
dB_ 

~dl ~ 



(9.18) 



g sm( 



mi p + m 2 a 
niip 2 + m 2 a 2 



Note that the last three equations are separated from the first four. This occurs in general 
when the Lagrangian active forces do not depend on the position of the point, but only on their 
velocities. For equal masses mi = m 2 , we have a further simplification: 



p sml 



dxi 

— = P cos9, 

dyx 
dt 
dx 2 

—— = a cos 9, 
dt 

dy 2 . 
— — = a sin 9, 
dt 



dp 

~dt 

da 
~dt 
dB_ 

dt 



g cos 9, 
■ g cos 9, 

- g sm0-5-; — 2- 
p z + a 



(9.19) 
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(ii) Second method for the single constraint equation (|9.14p . In this case, 

[Ci] = \m , - x 2 , - yi , xi] 

does not have the maximal rank for vi = V2 = 0. This is a singular state for whatever 
configuration (see Remark l7.3p ; the set of the singular states is the zero-section of TQ. Moreover, 
since 



we have: 
[C l 
G-- 

\ci 



£ 

o ^ o 

4 










1 

m 2 



mi 

G- 1 



+ 



m 

mi 
±f + yf 



m 2 

m x 2 

mi ' mi 



x 2 

mi ' 
2K 

mi m 2 ' 

_ ii 

m 2 



yi xi 

777-2 ' m 2 
-1 



G~ 



mi m 2 
2K ' 



xi 
m 2 



TT 



'■I] 



G~ 



yl 


x 2 y 2 


mm 


xim 


m\ 


m\ 


mi m 2 


mi m 2 


x 2 y 2 


x \ 


x 2 yi 


±i± 2 


2 

?77| 


2 

?7l| 


mi m 2 


mi m 2 


yim 


y\x 2 


yl 


xiyi 


mi 7772 


mi m 2 


2 

m 2 


2 

mi 


xiy 2 


Xl± 2 


xii/i 


x\ 


mi m 2 


mi m 2 


mj 


m\ 



7T" 



1 



















7771 











1 










mi 












1 









777 2 













1 








777 2 







^22/2 


mm 


xim 




2 

m\ 


2 

r?7| 


mi 7712 


mi 7772 




X 2 i)2 


x l 


X2ijl 


±i±2 


7771 Tn>2 


2 

m\ 


m\ 


mi 7772 


mi m 2 


2K 


mm 


mx 2 


iA 


xim 




mi 7772 


mi m 2 


ml 


ml 




x\y 2 


Xi±2 


xim 


x\ 




mi 7772 


mi 7772 


2 

mi 


2 

mi 



11 



-.12 



TT 



11 



1 



777 2 i) 2 



.12 



mi 2K mi 

mi 7772 ±22/2 



— 1 



2K 



m 



1 

7771 
7772 X 2 y 2 

2miK' 



m 2 y\ 



2K 



2K — m 2 yl 
2miK 



mi (xl + yf) + 777 2 xl 
2 mi K : 
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,13 



14 



,22 



,23 



,21 



,33 



,31 



7T 



13 _ mi m 2 ym 



VW2 



IT 



11 



7T 



22 



7T 



23 



7T 



21 



2K mi m 2 2A" ' 
_ mi m 2 xiy 2 _ ^12/2 
2if mi m 2 ~ 2K ' 
1 m 2 x\ 
mi mi 
_ mim 2 x 2 yi _ x 2 y\ 
2K mi m 2 ~ 2K ' 
mi m 2 x\x 2 x\x 2 



in 

mi 



m 2 x 2 
2K 



2K — m 2 x\ ra\ (xf + y\ ) + m 2 y 2 



2m : K 



2m 1 K 



2K 



7T 



7T 
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2X mi m 2 

J_ _ mi _2/f_ . 
m 2 2X m 2 
mi m 2 ±iyi m\Xiyi 



-Li 

m 2 



2K mf 



.9 



n 



7T 



11 



1 

m 2 



2m 2 ^ 
1 

2X m 2 m 2 



mi xf 



1 



"ii2/i 
2A" 



mi xf 
2K 



2ET — mi iji m 2 (x 2 , + y 2 ) + mi xf 



2m 2 K 



2m 2 K 



2K — m\ xf m 2 (x 2 + y 2 ) + m\ y( 



2m 2 K 



2m 2 K 



M = [0,0, 0,0], [d j cci] = [o\. 

As in this case Li = Ai, we have 

D i ■ {g ij _ ^ L . _ ( y djC C i = {g ij _ ^ L . = {g ij _ ^ y 



3i 



D 

D 2 

D 3 



1 
2K 

1 
2K 

1 
2K 

1 
2K 



A 



mi (if + ij 2 ) + m 2 x 2 



+ A 



m 2 x 2 y 2 



m 2 x 2 y 2 

Ax — + A 2 

mi 



-4iyiy 2 - A 2 ± 2 yi + A 3 



mi mi 
mi (xj + yf) +m 2 y\ 



mi 



+ A 3 yiy 2 - A 4 x x y 2 , 



A s x 2 y\ + A 4 ±ix 2 



m 2 (x 2 + y 2 ) + mi xf "iixiyi\ 



m 2 



A ■ ■ , A ■ ■ i A m l X l^l , A 

A±xiy 2 + ^4 2 xix 2 + A 3 h A± 



m 2 



+ A A 

m 2 J 

m 2 (x| +yl) +m 1 yj \ 

m 2 J 



Then the dynamical system (|8.12p reads 

' A x 



dx\ 
~dt~ 
dyi 

dt 
dx 2 

~dt 

dy 2 
dt 





dx± 


1 




~dT 


~ 2K 




dyi 


1 


= m, 


dt 


~ 2K 


= x 2 , 


dx 2 
IT 


1 

~ 2K 


= m, 


dm 


1 




dt 


~ 2K 



mi (xf + yf) + m 2 x 2 



mi 



+ A 2 
i _l .-,2 



m 2 x 2 y 2 
mi 



m 2 x 2 y 2 mi (xf + y{) + m 2 y^ 
A\ h A 2 



mi 



mi 



+ -43 2)12/2 - -44 x x y 2 ) , 

- A 3 x 2 yi + A A xix 2 

. . . m 2 (x| + yf) + mi xf miXiyA 

-4iyiy 2 - A 2 x 2 yi + A 3 h A 4 , 

m 2 m 2 / 

, ■ ■ , A ■ ■ , A . , m 2 {xl + yl) + m 1 yl 
- A\x\y 2 + A 2 xix 2 + A 3 h A 4 



m 2 



m 2 



For two points running on an inclined plane, 
mi (i;f + yf)+m 2 (xi + y 1 y 2 ) 



D 1 



D 2 



mi (xf + y\ ) + m 2 (x| + y|) 

wt2^2 (2)2 - m) 

"ii (if + yf) + ^2 (i| + y 2 ) ' 



L> 4 



m 2 (x^ + y|) + mi (xf + yijfg) 
mi (if + 2/f ) + m-2 (x\ + y\ ) 

"ii ^1 (2/1 - 2/2) 
mi (xf + yf ) + m 2 (x\ + yf) ' 
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For equal masses, mi = 777-2 




x\ + y\ + x\ + y\ ' 



±\ + yl + ±l + yl' 



Note 9.1. It is easy (and obvious) to propose examples of non-linear constraints: it is sufficient 
to choose any set of non-linear independent equations C a (q,q) = 0. However, any example 
of a non-linear constraint remains meaningless unless we know how to realize it physically by 
means of realizable devices. The famous Appell-Hamel example gives a matter of discussion 
(see [9], Ch. 4, § 2). Indeed, in order to be really a non- linear device, a certain distance of the 
Appell-Hamel device must be infinitesimally small. This fits with the thought of Hertz: non- 
linear constraints can be realized by passing to the limit x — > of certain physical quantities x 
(masses, lengths, etc.) in devices realizing linear constraints. 

The same kind of problem arises in trying to 'realize' two mass-points moving with parallel 
velocities. A tentative project has been presented in [I]. In fact, for an effective project, we 
have to invent devices for 

1. Realizing a mass-point. 

2. Realizing a parallel transport on the plane. 

3. Transforming the direction of the velocity of a point into a solid segment. 

4. Applying forces of special kind to the points (the weight is of course always present). 

This research is a work in progress. Updated information will be found on my personal web-page. 
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